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Abstract 

Third political parties are influential in shaping American politics. In this work we 
study the spread of a third party ideology in a voting population where we assume that 
party members/activists are more influential in recruiting new third party voters than 
non-member third party voters. The study uses an epidemiological metaphor to develop 
a theoretical model with nonlinear ordinary differential equations as applied to a case 
study, the Green Party. Considering long-term behavior, we identify three threshold 
parameters in our model that describe the different possible scenarios for the political 
party and its spread. We also apply the model to the study of the Green Party's growth 
using voting and registration data in six states and the District of Columbia to identify 
and explain trends over the past decade. Our system produces a backward bifurcation 
that helps identify conditions under which a sufficiently dedicated activist core can 
enable a third party to thrive, under conditions which would not normally allow it 
to arise. Our results explain the critical role activists play in sustaining grassroots 
movements under adverse conditions. 

MSC primary 37N99, 91F10; secondary 91D30 

Keywords grassroots, third parties, dynamical systems, backward bifurcation, threshold 
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1 Introduction 

The 2000 United States presidential election was for many a testimony to the impact of third 
parties in a traditionally bipartisan government. Ralph Nader, the presidential candidate for 
the Green Party, won 2% of the popular vote, a percentage that many attribute to the defeat 
of Democratic candidate Al Gore (Southwell, 2004). The Green Party captured a seemingly 
insignificant number of votes relative to majority percentages, yet its presence in the election 
ultimately served to shape American politics for the years following. This incident demon- 
strates how third parties, often emerging as grassroots movements (i.e., movements at the 
local level rather than at the center of major political activity), can ultimately impact at 
the national level, hence prompting the need to study their emergence and spread within a 
voting population. 

Third parties are defined as political parties operating along with two major parties 
in a bipartisan system over a limited period of time (where we define a limited period of 
time as a range of a few years). For the purposes of this paper we apply this definition 
to all minor parties. Traditionally, third parties have served as venues of political dissent 
for voting individuals dissatisfied with the major candidates in an election. They often 
tackle specific issues otherwise ignored by major political parties, thus relinquishing popular 
support nationwide. As Supreme Court Justice Earl Warren wrote in 1957, "History has 
amply proved the virtue of political activity by minority, dissident groups, which innumerable 
times have been in the vanguard of democratic thought and whose programs were ultimately 
accepted" (Warren, 1957). Hence, while third parties rarely capture the majority vote, their 
agendas, often incorporated into major party platforms, are significant nonetheless. 
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Given the potential relevance of third parties to national politics, we study, qualitatively 
and numerically, the dynamics of the emergence and spread of third parties on a local level 
where growth is measured in terms of the number of third party voters and members. We 
restrict our study to a local level because third parties usually originate in a small group and, 
via a "bottom-up" method of diffusion, spread within a population by acquiring local official 
positions and then expanding to higher levels of government (Kowalewski, 1995). Although 
individual personalities and circumstances dominate the initial formation of any group, the 
ability of even a small group of people to make itself heard within a larger community is 
great enough, especially in an age of information technology, that we consider the local level 
large enough to be described by a collective average. That is, the voting population under 
study is large enough that stochastic (random) effects are dominated by the deterministic 
average behavior of the group. 

We use an epidemiological paradigm (Brauer and Castillo-Chavez, 2001) to translate 
third party emergence from a political phenomenon to a mathematical one where we assume 
that third parties grow in a similar manner as epidemics in a population. We take this 
approach following in the steps of previous theoretical studies that model social issues via 
such methods (Bettencourt et al., 2006; Castillo-Chavez and Song, 2003; Crisosto et al., 
2001; Gonzalez et al., 2003; Song et al., 2006). The epidemiological metaphor is suggested 
by the assumption that individuals' decisions are influenced by the collective peer pressure 
generated by others' behavior; the "contacts" between these two groups' ideas are analogous 
to the contact processes that drive the spread of infectious diseases. Here we assume that 
a certain subpopulation of the voting individuals, defined according to certain demographic 
factors, is more receptive (in epidemiological terms, susceptible) to third party ideology than 
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the rest of the voting population, and that their political behavior is therefore driven by such 
collective peer-pressure contacts (Southwell, 2003; Timpone, 1998; Wong, 2000). 

There are many components that may affect the decision of a person when it comes to 
voting. During a political campaign all candidates spend a great amount of effort to assure 
that as many people as possible get exposed to their ideas, and most importantly, their 
name. People are often so overwhelmed by exposure to different candidates that they end up 
depending more on informal ways of obtaining information such as talking to other people. 
This makes networking among people very influential during election time; in fact some 
studies show that voting is "contagious" (Nickerson, 2008). It has been shown that people 
often rely on friends, relatives, coworkers, etc. to obtain information about the candidates 
(Robinson, 1976). Also, there is evidence that people who support a particular candidate 
tend to encourage others to vote for that candidate (Huckfeldt and Sprague, 1995). This 
makes the interaction between people a very important factor in the voting process because 
both the person wanting information seeks it from other people and the ones with particular 
preferences want to share them. 

Collective behaviors such as voting have been studied for decades, notably in the seminal 
work of Granovetter (1978), who considered individuals in a population to have a distribution 
of thresholds with regard to their willingness to participate in a particular collective behavior 
(he used rioting as a primary example). Other subsequent studies considered the role played 
by a core group of especially influential individuals (Macy, 1991; Oliver and Marwell, 1988), 
such as party activists recruiting voters. These studies used probabilistic and stochastic 
models which yielded primarily numerical results; in applying the epidemiological framework 
described above, we focus on the collective (rather than individual) thresholds for persistence 
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of a behavior (here voting and party membership) which our deterministic models allow us 
to calculate. 

While our model is designed to pertain to all third parties, we consider the Green Party 
as a case study. Although formally united under the Association of State Green Parties in 
1996 (and later the nationalized Green Party of the United States in 2001), state-based green 
parties have thrived in the U.S. at the local level since 1984, when the Green Committees 
of Correspondence (CoC) were formed with the purpose of organizing local Green groups 
and working toward the founding of a national Green political party (Marks, 1997). Our 
particular study focuses on the growth (and in some cases decline) of the Green Party in six 
states and the District of Columbia in the past decade, using voting and registration data. 
In comparing the predictions of our model to particular data, we consider a short time frame 
so that we can assume that social structure within the state in question does not change 
drastically, a necessary condition for assuming voting population heterogeneity. 

We organize our paper as follows: Section 2 describes our theoretical model, as well as a 
simplification in which the entire population under study is equally receptive to third party 
ideas. Section 3 presents the mathematical analysis of the simplified model (which employs 
qualitative analysis techniques from the field of nonlinear dynamical systems, as well as 
sensitivity analysis) and interprets the results, and Section 4 applies them to a case study, 
using data to estimate model parameters. Section 5 performs a similar but more limited 
analysis of the more complex model, and Section 6 draws conclusions about the implications 
of our models for the growth of third parties from grassroots movements. 
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2 A Population Model for the Spread of a Third Party 



2.1 Underlying Assumptions 

In developing our model, we apply epidemiological terminology to describe the growth of 
a third party. The assumptions we make about how individuals behave, and change their 
behavior, define the classes in our compartmental model and the rates at which individuals 
move between classes. 

(1) We assume that our population is a heterogeneous mix of individuals who belong to 
different backgrounds according to certain demographic factors. 

Our model considers a population of all voters, N, divided into two classes or sub- 
populations whose susceptibility to third party ideology is based on demographic factors such 
as education, socioeconomic status, race, gender, age, political orientation and professional 
occupation. Inherently, certain demographic characteristics, labelled as high affinity, make an 
individual more likely to subscribe to a third party's ideology, which targets a more specific 
audience than alternative majority agendas. That is, upon entering the voting system, 
certain individuals are more statistically inclined to vote a certain way. For example, a 
progressive environmental activist is statistically more likely to agree and vote for the Green 
Party agenda, which stresses communal-based economics, local government, and gender and 
racial equity, than a conservative corporate executive whose economic philosophy directly 
conflicts with that of the Green Party. For this reason we consider population heterogeneity 
vital to this study. 
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We apply the following method of dividing the entering voting population into two sus- 
ceptible classes: if an individual has more high affinity factors than low affinity factors then 
that person directly enters the high affinity class and similarly for low affinity susceptibles. 
We define high affinity factors as features of the individual based on his/her demographic 
profile that make him/her more inclined to vote for the third party; conversely, low affinity 
factors make the individual less statistically likely to subscribe to the party's platform. We 
assume that a constant proportion p (0 < p < 1) of new voters enters the high- affinity class 
H (the remaining 1 — p proportion enter L) . 

(2) We assume that individuals' affinity factors (demographic characteristics) remain fixed 
for the period of the study. 

(3) We assume that a third party's agenda remains consistent over time. 

We assume that individuals do not move from one susceptibility class to the other. One 
reason for this is the relative permanence of individuals' demographic characteristics. We 
limit our model to tracing the expansion of the third party; hence, we refer to a shorter time 
period over which we assume social structure remains constant. In other words, individuals 
with high affinity to the third party do not become individuals with low affinity, and vice 
versa. 

The other reason for this consistency has to do with the parties themselves. In addition 
to being more specific than major party agendas (i.e., more specific in their goals and less 
geared to moderacy), third party platforms tend to be more consistent over time. Third 
parties are not pressured to constantly adjust to the shifting demands of the populace since 
they do not seek the majority vote. Consequently, they do not target the majority voting 
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population. Each party has its own agenda, which appeals to certain sectors of the voting 
population. Hence different parties target voting populations that are more inclined to sub- 
scribe to their ideology. While one party, for example, may target individuals from a certain 
educational background that we, in our model, label as high affinity and that other parties 
may overlook, all parties nonetheless recognize that education factors into an individual's 
likelihood to support or refute that party's platform. It is true that individuals from varied 
backgrounds comprise the main parties, yet, when dealing with the specific agendas of third 
parties that do not strive to sway the majority vote, we assume that third parties appeal to 
individuals of certain demographic backgrounds more than others. Therefore, we account for 
the aforementioned standard set of demographic factors that parties look at when spreading 
their ideologies. In our paper we apply our model to an individual case study of the Green 
Party of Pennsylvania; however, the same methodology of distinguishing susceptibles can be 
applied to all third parties. 

(4) We define two levels of participation in third-party politics: voting for third-party 
candidates, and membership. A party exists only if it has members; we define members 
as those who pay dues, volunteer, and preside over party affairs. 

As described above, all voters enter the voting system either to the low affinity, L, or 
high affinity, H, susceptible class. According to our epidemiological metaphor, in addition to 
these two susceptible classes, our model includes three infected classes: Vh, Vl, and M, third 
party voters from the high affinity class, third party voters from the low affinity class, and 
party members respectively. We define party members as voters of the third party who pay 
dues to the party; often such members officiate, volunteer and actively campaign for voter 
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recruitment. In epidemiological terms, Vh and Vl correspond to voters of a lower degree of 
infection and individuals of the M class are voters infected to a higher degree. We distinguish 
between Vh and Vl because of their interactions with their respective "neighbors" in H and 
L. 

(5) We assume that third parties, emerging through resource-limited grassroots efforts, 
spread primarily via primary (direct) contacts between third-party supporters and 
susceptibles. 

(6) We assume that third party members have a greater effect upon voter recruitment than 
do third party voters, due to members' activism. 

In our model, a system of nonlinear differential equations, we consider susceptible move- 
ment into voting and member compartments as well as possible regressions back from the 
third party voting phase into the susceptible class. Once an individual is susceptible he/she 
can become 'infected' (either Vh or Vl) through direct contact with the Vh, Vl, and M 
classes. Due to a lack of funding and resulting lack of mass-media exposure to the gen- 
eral population, these primary contacts with susceptibles involve such direct interaction as 
personal meetings, phone conversations, and electronic communications like personally ad- 
dressed e-mails and weblog (blog) comments. We assume that the rate at which these con- 
tacts occur is proportional both to the size of the susceptible group and to the proportion 
of third-party supporters in the population. 

We do not consider a linear term weighing the influence of media coverage from the third 
party (i.e., secondary contact factors) in the forward transition from both susceptible classes 
to third party voting classes. Instead, we focus on the nonlinear terms considering the effects 
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of voters from the Vh, Vl, and M classes, where voters from Vh and Vj, bear an affinity- 
specific influence (3h (from H to Vh) and (3l (from L to Vl) in third party voter recruitment. 
Through activism, members from M influence susceptibles of each type at higher rates {af3n 
and a/?L, respectively) than voters, with their increased influence measured by the parameter 
a (a > 1). 

(7) We consider both primary (direct) contacts and secondary (indirect) contacts in the 
regression of third party voters to the susceptible class. 

(8) We assume that all other parties exert equal influence in discouraging third party 
voting. 

(9) We assume that individuals have more influence upon others of their same affinity 
class (high or low) than upon members of the other affinity class, in encouraging and 
discouraging third party voting. 

We consider the transitions back from the third party voting to the susceptible classes to 
involve both linear terms, e#V# and €lVl, and nonlinear terms, 4>h(H + uL)^- and 4>l{L + 
aH)j^, contributions by secondary contacts with the opposition (i.e., media from well-funded 
majority voters) and direct contact with the susceptible classes respectively. Compared to 
primary contacts, described in a previous paragraph, secondary contacts include mass e- 
mails, media, and circulating literature. 

Voters from a certain susceptibility class (with its own set of demographic factors) address 
issues that usually appeal more to voters deriving from the same class. Therefore, susceptible 
individuals with high affinity bear a greater influence in recruiting voters who came from 
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the high affinity susceptible class back into the susceptible class than susceptible individuals 
with low affinity. The reduction in influence by individuals from a different affinity class is 
denoted by a, so that in regressing back from Vh to H, oL represents the lesser influence that 
L individuals exert on voters from Vh than do susceptibles from H, the higher affinity class 
(similar reasoning applies to the Vi-to-L transition). Likewise, the cross-affinity influence in 
recruiting third party voters is reduced by a factor of er. 

(10) We assume that third party voters become active party members through the ongoing 
efforts (primary contacts) of the members, who have made a permanent commitment 
to the party. 

Once voting for the third party, individuals can become party members. They enter this 
higher state of infection via the nonlinear terms jVh^ an d 'jVlj^, where we only consider 
the influence of primary contacts with party members in bringing about this transition, 
measured by the rate parameter 7. Given that we are studying the spread of the party, 
we assume that party members do not resign their memberships. We reason that once an 
individual feels strongly enough to join a party, he/she retains his/her loyalty to the party; 
the only way a person stops being a member (during the growth period under study) is by 
leaving the voting system. 

Finally, we consider natural exits from all classes as a result of death or moving. The sum 
of the equations of the model, for both versions developed below, gives ^ = 0, verifying that 
our population stays constant, a safe assumption by which the number of people entering the 
voting system (i.e., coming of age, moving in) counterbalances the number of people leaving 
the system (i.e., dying, moving out). 
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2.2 The General (Two- Track) Model 

Ph(Vh + <tVl + oM)/N 




(l-p)fiN 



v„ 


\ 

V 










Figure 1: The two-track model, with per capita flow rates 



We first introduce a two-track model, as described immediately above, to study the 
dynamics between a heterogeneously mixed population of susceptible voters, third party 
voters, and party members. We apply the following set of ordinary differential equations to 
model voting dynamics, as illustrated in Figure [TJ 



^ = pfjiN + e H V H + <f) H (H + crL)jj-- (3 H (V H + aV L + aM)^ - fiH, (1) 

^ = (1 - p)fiN + e L V L + ML + *H)^- (3 L (aV H + V L + aM)^ - /iL, (2) 
at N N 

^jf = (3h(V h + aV L + aM)-- e H V H - <f> H {H + <tL)^§ - ^-p- - yV m (3) 

at JM N Jy 

-j± = f3 L (aV H + V L + aM)- - e L V L - L (L + aH) V -± - - ^V L , (4) 

N = H + L + V H + V L + M. (6) 



Adding equations (1), (2), (3), (4) and (5) yields ^ = 0, showing that the total population 
TV is constant over time. Model parameters are summarized in Table 1. 
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2.3 The Simplified (One- Track) Model 



In order to facilitate analysis of the two-track model, we initially consider a simplified version 
that does away with voting population heterogeneity (assumption (1), and consequently (9)) 
before exploring analysis for the more complex system. This simplified model assumes a 
homogeneous susceptible population (p = or p = 1), reducing the two-track model to one 
susceptible class, S, and two infected classes: third party voters, V, and party members, M, 
respectively. The S class comprises those individuals who vote, but do not vote for the third 
party. The V class comprises the third party voters, and the M class again has third party 
members (i.e., party officials, donors, volunteers). 

In the one-track model we omit unnecessary parameters from the heterogeneous version. 
Figure [2] illustrates the one-track model, and Table [2] summarizes the parameters. 



3(V + aM)/N 




7M/N 




Figure 2: The one-track model, with per capita flow rates 

In this case the model reduces to the following system, which is effectively two-dimensional 
since N can again be seen to be constant: 



— = fiN + eV + <PS- - f3{V + aM)-- »S, 



V !MV 
eV-cj>S-- nr - f iV, 



(IM ^MV 



/iM, 



dt N 
N =S+V+M 



(7) 
(8) 
(9) 
(10) 
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3 Analysis of the One- Track Model 

We begin our analysis by calculating equilibria for our model and determining conditions 
for their existence and stability. We first simplify the system in two ways. Since the total 
voting population iV is assumed constant, we can reduce our system to two dimensions by 
rewriting S = N — V — Mm equation (jSJ), so that 

dV N -V - M V >yMV 

^ = (3{V + aM) V - eV-^N-V-M)^-^--^. (11) 



We can now analyze the system defined by ffTTj) and ([9]), as S can always be found once V 
and M are known. 

The constancy of N also allows us to proportionalize the system by defining new variables 
v = m = and s = = 1 — v — m, which give the proportion of the population in 
each class. Dividing equations (ITTj) and ([2D by N and substituting the new variables gives, 
finally, the system 

dv 

— = f3(v + am)(l — v — m) — ev — 0(1 — v — m)v — ^mv — fxv, (12) 

(JjL 

In order to analyze stability we linearize the system and compute partial first derivatives 
with respect to each of the variables, v and m, obtaining the Jacobian matrix J\ for system 

flU-GHD: 



Ji 



(/3 - 0)(1 - 2v - m) - (a/3 + 7)m - (// + e) ap(l -v- 2m) + 7)^^ 



7m — 
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3.1 Ei: Party- Free Equilibrium (PFE) 

The party-free equilibrium (PFE) for the reduced system occurs at (0,0), the steady state 
achieved when the entire population resides in the S class (i.e., the third party has neither 
voters nor members and, by definition of party existence, does not exist). The PFE is 
essentially analogous to the disease-free equilibrium in epidemiology and always exists as a 
possible outcome for the voting population. 

Applying the above reduced Jacobian matrix to our PFE, (0,0), where we only consider 

the v and m terms, we determine PFE stability: 

/ \ 

{f3 - 4>) - (// + e) af3 

Ji(0,0) = 

V -"J 

The equilibrium point (0,0) will be locally asymptotically stable (LAS) if all the eigenvalues 
of the matrix are negative. Assuming jj, > 0, the eigenvalue — /i of the Jacobian is always 
negative, whereas the second eigenvalue {(3 — 0) — (fi + e) < if and only if < 1. For 
ease of notation and interpretation we define the threshold quantity R\ = , so that the 
PFE is LAS if and only if Ri < 1. We discuss the relevance of this threshold value in the 
last part of this section. 

3.2 E 2 : Member-Free Equilibrium (MFE) 

The member-free equilibrium (MFE) occurs when M = but V, S ^ 0, i.e., the voting popu- 
lation subdivides between susceptibles, S, and third party voters, V. While mathematically 
possible, this outcome is politically unrealistic given that voters cannot vote for a party 
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that does not exist (recall our assumption (4) that party existence depends on the presence 
of an M class). For mathematical consistency, however, we consider the equilibrium point 
(v^m^), where = 0. (Note that the asterisk superscript denotes equilibrium values, while 
the numerical subscript distinguishes the equilibrium point in question.) This arises from 
the equilibrium condition obtained by setting dm/dt = in ( fl3l) : 

(7^* - n)m* = 0, (14) 

which implies that either {^v* — fi) = (which we will consider later) or m* = 0. 
To find i>2 we set dv/dt = and m* = in ( TT2|) and rearrange terms to get 

(J3 - 0K 2 + (// + t + - (3)v* = 0. 

This implies that either v* — (the party- free equilibrium) or (j3 — <ft)v* = (3 — (fi + e + (f>). 
We consider the situation where v* 7^ 0, solve for v* and simplify the results as follows: 

Vn = 1 = 1 

2 ~ (f> R X 

where v% retains political value only if Ri > 1 — otherwise v % < which is meaningless. 

Finally, we can also write = 1 — v% — m\ = 1 0^ > = which makes sense politically 
only if Ri > 1 (i.e, < 1). We can therefore express our member- free equilibrium as 
E2 = (1 — 5^, 0). The MFE exists if and only if R 1 > 1, since ignoring this condition leads 
to a negative third party voting population. 

The above situation makes mathematical sense but not political sense since parties, by 
our original assumption, do not exist without members, and in this member-free case we deal 
with voters who vote for a non-existent party. We might interpret this situation as having 
voters still willing to vote for this party, but no party candidates for whom to vote. 
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Regardless of the political likelihood of MFE existence, we consider its stability. Again, 
we apply the method of using the reduced system's Jacobian matrix in determining the 
stability of the member-free equilibrium: 



■Mi-ib°) 



(ji + e)(l - Rt) a(3(jr) - (/3 - + 7 )(1 - 



V 7(1 A* / 



The reduced system equilibrium point (1 — ^-,0) is locally asymptotically stable if all the 
eigenvalues of the above matrix are negative. We know that, since R\ > 1 (in order for the 
MFE to exist), one of the eigenvalues, (p+ e)(l — Ri), is negative. The second eigenvalue of 
the Jacobian is 7(1 — -^-) — //. This eigenvalue is negative if and only if —(1 — -^-) < 1. 

We define the left hand side of the inequality as R2 = ^(1 — pz$) — -(1 ~ jjr)- Hence 
we have derived two threshold parameters Ri and R 2 that determine equilibria stability 
depending on relative parameter values. (Note that they are related: Ri = 1 4^ R% = 0.) 

3.3 E% and E41 Survival Equilibria 

In the event of survival equilibria, the voting population subdivides between susceptibles, S, 
third party voters, V, and members, M. We regard this as a successful state of coexistence 
and, given certain conditions, the point at which the party thrives. We determine the 
equilibrium proportions by returning to the condition (|14p that dm/dt = 0. Since ^ = 
(jv* — /i)m*, v* = ^ when m* 7^ 0. Here we impose the condition \i < 7 so that v* < 1, 
since s* + v * + m* = 1. 
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Next we set dv/dt = and v* = in (pj: 

^. = (/3- < f,)(L (i-^- m *\ + a p m * ( i _ t - m * ) _ + e + 7m *) ^ = o. 
di 7 V 7 / V 7 / 7 

This equation, which is quadratic in m*, can be rewritten (after dividing through by —a/3) 
in the form /(to*) = m* 2 + £?M* + C = 0, where 



B 



(3- <p\ fi 



a(3 I 7 



+ 



a/3' 



C=- 







a/3 



1 - 



/i\ yU + e 



7 



a/3 



(15) 



Solutions to /(to*) = are given by the quadratic formula 



m i 



2 L 



however, in general, these solutions may be real, or may fall outside the meaningful interval 
(0, 1 — v*\. Depending on the values of model parameters, there may be 0, 1, or 2 solutions 
within this interval, each corresponding to a meaningful equilibrium with positive party 
membership. 

Analysis of the conditions involved in determining the number of survival equilibria is 
considerably more involved than that for E\ and E 2 ] details are given in the Appendix. The 
results, which introduce an additional threshold quantity, can be summarized as follows. 



Proposition 1. (i) If R 2 > 1, the system (fl2l -( TT3|) has precisely one survival equilibrium 

(ii) If R 2 < 1, then the system has two survival equilibria, E 3 and E4 = (^,m*_), if and 
only if R3 > 1, where 



r-2 



R 3 = mm(R 3a ,R 3b ), R 3a = r 3 ( 1 - ^ ) , R 3b = ^ ( 1 - J^—^- + h ) , (16) 



r 2 
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and 



_ r ^ = a pi h _ 2 | A + £ _ 1 

a/3 v 7 ^ V V 7 



Otherwise there are none. 



To interpret the conditions i?3 a > 1 and i? 3 {, > 1, we can rewrite them as follows: 



B < & R 3a > 1 «»— + ^ < - , (17) 
ap 7 /i 



B 2 - AC > & R 3b > 1 ^ \ + J ^—2 + fc(e/ T ) < . /- , (U 

V a/3 V 7 V ^ 



where 



*<^=V^vKv r ^" i ) >0 - 

Both inequalities (|T71) and (|T8|) relate the minimum amounts of time taken for an indi- 
vidual to be influenced by a party member (M) to move from S to V, l/a/3, and from V to 
M, I/7, to the average lifetime of an individual in the voting system, In order for the 
party to survive, a weighted sum of the first two times must be less than the average lifetime 
in the system, or else individuals in S will, on average, leave the system before they can 
"replace" the members in M who recruited them. The weights in the sum involve q = 
the relative effectiveness of voters V to members M in recruiting new voters from S, because 



the influence o 
to some extent 



jparty voters reduces the recruiting threshold burden on party members M 
i (Note that < q < 1.) 



Observe that (IT51) implies that 



a/3 V 7 V A* 



1 If we define 7\ = l/a/3, T 2 = 1/7, and T 3 = l//i, (TIT)) and ([T5]) can be written more simply as 
Ti + (1 + q)T 2 < T 3 , VT\ + ^{l-q)T 2 + h < y/Ti. 

20 



If e << 7, we can generate a Maclaurin (Taylor) expansion in e/7 for (|T8|) which yields 

2" 




a/3 + V 7 + 2 



— /- 

a /3 V A 4 



7 



< 



This expansion suggests more quantitatively how the influence e of voters for other parties S 
in getting third-party voters V to "defect" back to the major parties complicates the third 
party's survival (which would otherwise only require (fT9l) ). 

Having established conditions for the existence of the two survival equilibria, E% and £4, 
we examine their stability. The reduced Jacobian matrix for the survival equilibria follows: 



(/?-</>) 1-2 



m. 



-(a/3 + 7)774 - (/i + e 



a/3 (l - H _ 2m* ± 



y 7"4 y 

The stability (LAS) criterion that the eigenvalues of this matrix have negative real part is 

equivalent to the conditions that the trace be negative and the determinant positive. We 



calculate 



det Ji{^m* ± ) = -7«4 



a/3 1 1-^-2774) - (R -11 

7 / 7 



so that 



det Jx(^ mi) > & 



a/3 I l-^-2mI ] - (B - <f>)^ - a 
1 J 7 



< 



7 



/3 — 0\ /X /i 



a/? / 7 a/? 



5 
"2 " 



Since m*_ < — < m* + , we see that £4 is never stable, while the stability condition for £3 
reduces to tr Ji(-,m+) < 0. Thus we calculate 
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tr Jtf, m* + ) = (/? - 0) ( 1 - 2^ ) - (ja + e) - (J3 - + af3 + j)m* 



+ • 
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so that 

(J3 - 0) (l - 2*) - (// + e) (/3 - 0) m (i? 2 - 2 ) 
7 (p - <p + ap + 7) (p — <p) + «p + 7 

This is true for all m* when R2 < 2. The case when R2 > 2 requires further algebra and 
is relegated to the Appendix. The result is that E3 is always stable (LAS) when it exists, 
while £4 is always unstable. 

3.4 Global Behavior 

Table [3] summarizes the conditions for existence and local stability of the four equilibria of the 
one-track model. Because system dZJ)-® can be reduced to a set of two differential equations 
( I12p - (|13p . we can apply the Poincare-Bendixson Theorem to establish global stability. Since 
the system is well-posed, with the invariant set D = {(v , m) : v, m > 0; v + m < 1}, there 
are no unbounded solutions beginning in the state space, and a straightforward application 
of Dulac's Criterion with the function b = 1/vm confirms that there are no limit cycles, 
either: 

d (^dv\ d f^ m \ < g • d 
dv \ dt J dm \ dt ) 

Therefore, all solutions to the system which begin within the state space must approach 
an equilibrium. In particular, when only one locally stable (LAS) equilibrium exists, that 
equilibrium must in fact be globally stable (GAS). 

The conditions in Table [3] can be graphed to show the different possible global behaviors 
of the one-track model. Each condition corresponds to a curve dividing parameter space into 
multiple regions, each of which represents a different global behavior (see Appendix IA.1I for 
the derivations). The five resulting regions are illustrated in Figure [3] and summarized in 
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Figure 3: Regions of equilibrium stability, with q, x, y as in Appendix A 

I In region I, E\ is the only stable equilibrium; the party will always go extinct. 

II In region II where R\ > 1, R 2 < 1, and R 3 < 1, E 2 , the member-free state, is the only 
stable equilibrium. 

III In region III where both Ri > 1 and R 2 > 1, E% is the only stable equilibrium; the 
party will inevitably approach a survival state. 

IV In region IV where R\ > 1, R 2 < 1 and R3 > 1, E 2 and E% coexist as stable equilibria 
although, if placed in a political context, E 2 is not a realistic outcome for the party. 

V In region V where R\ < 1, R 2 < 1, and R3 > 1, E\ and £"3 are both stable; depending 
on the initial conditions the solution tends to one state or the other. 

The coexistence of two locally stable equilibria (one representing party survival and the 
other party extinction) in regions IV and V involves a phenomenon known in epidemic 
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modeling as a backward bifurcation, in which the exchange of stability at a threshold value 
(here, R2 = 1) reverses direction, creating a situation in which the unstable equilibrium E± 
serves to split the state space into two basins of attraction for the two stable equilibria. In 
these regions, the equilibrium approached depends upon initial conditions: in particular, a 
large enough core of dedicated members M can sustain the party (toward E3). What is 
unusual here is that the "backward" part of the bifurcation curve can extend back beyond 
not only the bifurcation at R% = 1 but also the bifurcation at Ri = 1 (where R2 = 0), as 
illustrated in Figure HI Whenever Ri < 1, R 2 < 0, a condition that would normally lead to 
the death of the party given local asymptotic stability of the PFE, there are still conditions 
(i?3 > 1) under which two survival equilibria exist, one of which is stable (E$). In other 
words, the party can thrive in conditions under which it would normally die out, given that 
we have the necessary parameters and sufficient initial number of M individuals. 

The ideal conditions for the party, of course, are when R 2 > 1 (region III), in which the 
party survives regardless of initial conditions. (Note R2 > 1 implies R± > 1.) The last part 
of this section interprets all these mathematical thresholds in political terms. 

3.5 Threshold Parameters Ri, i? 2 , and R 3 

Our system contains three local thresholds or tipping points where population outcomes, 
measured as S, V, and M, depend on parameter values. By tipping point we refer to the 
sociological term that describes the point at which a stable phenomenon turns into a crisis, 
which, in a political context, corresponds to the extreme states of the party: death and 
growth (Gladwell, 2000). In the context of our model, for example, the party cannot sustain 
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Figure 4: Diagram illustrating a backward bifurcation that extends below R 2 = (Ri = 1), 
using parameter values = 0.15, \x = 0.05, e = 0.25, 7 = 0.20 (all rates in yr" 1 ), a = 1.25 

grassroots growth until parameter conditions reach R± — 1, after which point the third party 
voting and member classes gain individuals. These threshold quantities correspond to similar 
terms in demographic and epidemic models called reproductive numbers which measure the 
average number of offspring or infections caused by a single group member during its lifetime. 
In our political context it is more appropriate to interpret these reproductive numbers in 
collective terms as regards the influence of party voters on the susceptible population. We 
distinguish between the aforementioned thresholds, R±, R2, and R3, by analyzing them 
qualitatively in a political context. 

Ri denotes the average number of susceptibles influenced to vote for the party by a single 
party voter in V, if dropped into a homogeneous population of susceptibles. The expression 
Ri = gives the net peer pressure, (3 — (f), on susceptibles S by voters V, multiplied by the 
average time, spent in the voting class V. As mentioned above, it is most appropriate 
to interpret Ri as an average number per voter taking into account the collective influence 
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of all the party voters. Note that we are assuming that the influence of party voters on 
susceptibles (3 is stronger than the reverse influence of susceptibles to discourage party 
voters, a necessary condition for party emergence, which guarantees that R± > 0. 

R 2 measures instead the second stage of recruitment: how effective party members M 
are in recruiting third party voters to become members once there are enough individuals 
in V. The expression R 2 — ^(1 — ^-) gives the product of the rate at which party activists 
M recruit voters from V into M, 7, the average political lifetime of a party member, -, 
and the proportion of the population N in V at the MFE, (1 — j^) (since only voters V 
can be recruited directly into M). Similar to R±, R 2 measures the average number of V 
to M conversions per individual in the M class. If party voters are ineffective at influenc- 
ing susceptibles to vote for the party and R± < 1, then the value of R 2 will be negative; 
rather than interpreting this as a negative ability of party members to recruit voters into 
activism, however, it should be seen as an indication that the pool of party voters available 
for recruitment into party membership is not there, preventing normal party growth. 

i? 3 measures the extent to which party members M actively recruit susceptible individuals 
S into the voting class V. This activism, which sidesteps the traditional hierarchical structure 
of a party, is a key characteristic of growing grassroots movements, which often lack a political 
environment favorable to their growth in the traditional way outlined above (R 2 < 1). 
Because of the two conditions required mathematically for party survival when R 2 < 1, R 3 
is defined as the greater of two quantities, both of which involve party members' abilities to 
recruit voters, r 2 , and members, r 3 , over their political lifetimes, as well as the relative efficacy 
of voter peer-pressure influence to members' influence through activism on the susceptible, 
q. The general form ^(1 — l/r 2 ) present in both components of R3 parallels the form of 
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i?2 = ^(1 — 1/Ri), with members' potential r 2 for converting third-party voters into members 
reducing the constraint on their ability r% to produce those voters in the first place. Note 
that r 2 must be great enough in order for R 3 even to be positive: that is, only if the members 
can recruit other members well can their recruitment of voters sustain the party. 
We now list possible outcomes involving these thresholds and their implications: 

1. When Ri < 1 (regions I and V in Figure [3] and Table H]) the net influence of party 
voters V upon susceptibles is weak enough that grassroots emergence of the party is 
not possible. The hierarchical structure of party involvement precludes normal growth 
of party membership M when party voters are unable to replenish their own ranks 

< 1 implies R2 < 0). However, exceptionally, if the recruiting ability of party 
members M is great enough (R 3 > 1, region V), a sufficiently large dedicated core 
can sustain the party through its own efforts, despite the relative inefficiency or lack 
of influence of those merely voting for the party. This outcome illustrates the key role 
activists play in party survival during periods of adversity. 

2. When R\ > 1 each individual in V is converting, on average, more than one person in 
S into V, thus allowing the voting class V to thrive. If, in addition, R 2 < 1 (regions II 
and IV), we have a situation in which party activists are then unable to recruit from 
the voting class effectively enough to maintain the party core M (perhaps because the 
voting class is too small). Normally this would lead to an outcome in which a group 
disposed to vote for the party remains (in V) but the party core itself dwindles away, 
leading effectively to party extinction (with the ideologically closest main political 
party perhaps adjusting a platform to capture these votes). However, as with the case 
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when R\ < 1, a sufficiently large initial core can ensure the party's survival when 
the core is effective at influencing individuals to begin voting for the party (R 3 > 1, 
region IV). 

3. The condition R 2 > 1 explains the case where the party (V and M classes) grows 
normally, by recruiting members from the S and V populations, respectively Here 
the party voters are influential enough in garnering new voters from S (Ri > 1) that 
the party's survival does not depend on party activists' ability to recruit new voters 
directly (i.e., R 3 does not come into play). In epidemiological terms, this corresponds 
to a successful invasion: conditions are so favorable for the development of the party 
that it will become established even with a small initial group of members. 

In general, the survival of the party is determined by the interplay among the three 
recruitment processes involved in the model: from S to V by members of V, as measured 
by Ri, from V to M by members of M, as measured by r 2 (R 2 incorporates both of these 
first two processes); and from S to V by members of M, as measured by r 3 (R 3 incorporates 
all three processes). Party survival requires either that the first two processes be effective 
enough (R 1 > 1 and R 2 > 1) for a small grassroots effort to take hold, or else that an 
initial membership core be large enough, and the second and third processes effective enough 
(R 3 > 1), that direct recruitment by party activists can sustain its membership. 
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4 A Case Study: The Green Party 



4.1 Methods 

As an application of the simplified one-track model analyzed in Section [3X we used Green 
Party registration and voting records for six states (CA, ME, MD, NY, OR, PA) and the 
District of Columbia (DC) to study the growth of the Green Party during the past decade 
or so. 

We began by establishing basic demographic information for the target (study) popu- 
lation. A report from the Pew Research Center (2007) provided annual statistics on per- 
centages of Americans identifying themselves as confirmed Democrats or Republicans, lean- 
ing Democrat or Republican, or confirmed or leaning toward third parties, for the period 
1987-2007. We took an average of these percentages during the period 2000-2007, during 
which time approximately 25% of the voting population identified itself as leaning Democrat 
or Republican, and an additional about 11% as confirmed or leaning toward third parties 
(including the Green Party). We estimate our target population — those capable of being in- 
fluenced to vote Green — as all of the former group, and about half of the latter group (since 
the latter group also includes confirmed Greens). Thus we estimate that the target popula- 
tion consists of about 25%+6%=31% of the total voting population of each state. For each 
state, we averaged the voting population size over the time period of interest (which varied 
slightly from state to state, as detailed below) from voting records. We then normalized the 
voting and registration records as proportions of the target population in each case. 

To determine the replacement (or mortality) rate /i, we used the average 2003 life ex- 
pectancy at birth in the U.S. of 77.5 years given in Shrestha (2006), and the minimum 
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voting age of 18 years, to derive an average voting lifetime of 59.5 years, for an estimate of 
jj, = 1/59.5 yr ~ 4.58 x 1CT 5 /day. (A United Nations (2008) report gives an average U.S. 
life expectancy at birth of 78.3 years for the period 2000-2005, which yields a comparable 
estimate of /i = 4.77 x 10~ 5 / day.) 

The data used for this case study came from official voter registration records and elec- 
tion results from each state (see bibliography for sources). Voter registration records showing 
Green Party registration totals, given in some states as often as monthly, were fit to the size 
of the member class M over time. Since the dates for which Green Party registration data 
were given varied from state to state, the initial and final times did as well, but covered 
approximately the decade 1999-2008. Since voting data were given less often (in general in 
November of even-numbered years) these data were used, where available, as initial condi- 
tions U(0) + M(0). In other cases V(0) was estimated along with other model parameters 
as discussed below. Table gives the time periods modeled as well as the initial conditions 
used. The size of each class is given as a percentage of the total target population (the 
size of which is also given in the table). The rescaled Green Party registration data is also 
shown in Figure As can be seen in the graph, some states saw a noticeable change in the 
Green Party's trajectory following the November 2004 election, and so the data for these 
states (DC, CA, OR, NY) was broken into two subseries (three in the case of DC, which also 
underwent a visible change following the 2006 election), with different parameter estimates 
for each subseries. 

Data fitting used the program Berkeley Madonna 8.3 to obtain a least-squares fit to the 
Green Party registration data M(t) for the model parameters (3 — $, a/3, e, and 7, and, where 
necessary, initial conditions. {(3 — (f> and a [3 were estimated directly, rather than a, (3 and 
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Figure 5: Green Party registration data for the seven states (and district) used in the case 
study, shown as percentages of the target population. 

separately, in order to reduce the number of free parameters.) The estimation process was 
iterated until the estimates were not on target interval boundaries, and this optimization 
was carried out 10 times for each state (with different initial guesses) to obtain the best 
possible fit. The resulting estimates are also given in Table [51 These estimates were then 
used to determine the state of the system in each state by calculating reproduction numbers 
and predicted end states (equilibria); results are shown in Table [6j 



4.2 Results 

As seen in Figure EJ the data indicate a clear period of growth for the Green Party in the 
selected states during the first five years of the twenty-first century. However, for some states 
(DC, CA, OR, NY) the data also show a visible decrease beginning at the end of 2004. The 
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data indicate a clear change in the political landscape in these states following the November 
2004 election, in which the positive trends the Green Party had been seeing in recent years 
reversed, and voters left the Green Party, probably for the Democratic Party which saw gains 
in the 2006 and 2008 elections. The DC Greens rallied in 2005 but experienced this same 
trend reversal following the 2006 election. In each case, Table |5] suggests that the primary 
reason for this change was a significant jump in the attention the target population paid to 
the media and the activities of the two primary political parties, as evidenced by a marked 
(orders of magnitude) increase in the value of the parameter e (despite, in most cases, a 
simultaneous increase in the net peer voter influence (3 — <ft). 

The other three states in this study (ME, MD, PA), on the other hand, have continued 
to see healthy sustained growth in the Green Party over the past decade, although the 
growth in Maine appears nearly meteoric compared to the slow building up in Maryland 
and Pennsylvania. This difference is reflected in the much higher value of R\ for Maine (cf. 
Table [6]). The parameter estimates for Maine are much lower than for Maryland (cf. Table [5]), 
likely reflecting the lower overall person-to-person contact rate (Maine has no large urban 
areas, whereas Maryland's population is concentrated in them) as well as perhaps lower 
sensitivity to others' opinions, but the ratios reflected in the threshold quantities provide 
fertile ground for the Green Party's growth. Our one-track model predicts (Table E]) a slightly 
higher equilibrium level of proportional participation in MD than in ME, but it will take 
much longer to reach that end state (during which time political influences may change). 

Pennsylvania, meanwhile, appears to be in a more unusual situation, as the parameter 
estimates obtained from the data suggest the presence of a backward bifurcation as discussed 
in Section [3j Here the tipping-point thresholds (i?i, R 2 ) do not appear to favor the long-term 
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survival of the Green Party, but an initially large influx of Green voters would allow the party 
to thrive, even enjoying majority status among the target voting population. The reason for 
this phenomenon (also seen in the initial periods for DC) can be observed in the parameter 
estimates, which show a low general sensitivity (3 — <p to the passive influence of community 
peers, but a high sensitivity a(3 to grassroots activists (relative to other PA parameters), 
so that a sufficiently high core of party activists could sustain the party. Despite its initial 
similarity to Maryland in terms of the proportional data illustrated in Figure [5j the growth 
of the Pennsylvania Greens is indicated by the model to be part of a transient response — it is 
just that, in this system, transient responses can last on the order of decades, during which 
political influences can change significantly. 



5 Analysis of the Two- Track Model 

Having thoroughly examined and interpreted the one-track model, we now perform analysis 
on the original, heterogeneous two-track model. Since the added complexity of this model 
precludes a complete qualitative analysis, we shall make use of numerical analysis when 
necessary to demonstrate behavior analogous to that exhibited by the simpler model. 

First, we proportionalize the two-track model in the same way we did for the one-track 
model to get: 
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-77 = W + ChVh + <j>H(h + crl)v H - Ph(vh + crv L + am)h - fih, (21) 
at 

^ — (1 - + e L v L + + a7iW - /3 L (av H + v L + am)l - fil, (22) 
at 

= Ph(vh + crv L + am)h - e H v H - <j) H (h + al)v H - ^mv H - /iv H , (23) 

^ = f3 L (av H + v L + am)l - e L v L - <f) L (l + ah)v L - ^mv L - (j,v L , (24) 
dm 

—— = ^ttivh + "ymvL — ftm. (25) 

C£6 



5.1 The Party-Free Equilibrium (PFE) and R[ 

The two-track model, like the simpler version, includes a party-free equilibrium. Observing 
that v* H — v* L — m* — satisfies dvu/dt = dvi/dt = dm/dt = 0, we substitute into 
dh/dt = dl/dt = to find the PFE; E 1 is (p, 1 — p, 0, 0, 0). The stability of E 1 is again tied 
to the first threshold quantity, which we shall denote R[ (we shall use prime superscripts to 
denote thresholds for the two-track model) and calculate using the next-generation operator 
method (e.g., Castillo-Chavez, Feng and Huang, 2002), where v H) v L and m are considered 
the infective classes: 

R[ = 7, r HH + f'LL + V 7 ( r HH -r LL ) 2 + A r LH r HL ] , (26) 



where 



Thh = V — ; — 7—— and r LL = (1 - p) 



fjL + e H + (l- p)<y</)H ' f^ + e L + pa<f) L ' 

respectively, are the high-affinity and low-affinity analogues of R\ (that is, the influence of 

vh on h, and of vl on /), and 

r LH = p — ■ ■ — and r H L = (1 - p) 



/i + e L + pa<f) L '' ' fi + e H + (l- p)v<Ph 
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measure the cross-affinity influences (of Vl on h and vh on I). See Appendix B for calcula- 
tions. 

R\ can be interpreted as the average number of individuals converted into a third-party 
voter (in either voting class) by an individual in vh or introduced into a population where 
no one yet votes for the given party. The average is somewhat complicated, as it involves 
four different contributing influences, each represented by one of the four r's in (|26|) . Each 
of these component numbers has the same form as R\ for the one-track model (q.v.), but 
measures of conversion of high-affinity voters are multiplied by a factor of p, the proportion 
of the population which has a high affinity for the given party, while measures of conversion 
of low-affinity voters are multiplied by the proportion 1 — p of low-affinity individuals. In 
addition, cross-affinity influences are reduced by the factor a. Note that in the extreme cases 
p = and p = 1 R[ reduces to R\ . 

We can further interpret the expression for R[ in political terms by observing that 

m&x(r HH , r LL ) < R[ < m&x(r HH , r LL ) + y/r LH r HL (27) 

(again see Appendix B for details). That is, R[ is at least as great as each of the within- 
track conversion efficiencies, and exceeds the maximum of the two (presumably thh) by 
less than the contribution of cross-track influences. This latter contribution is the geometric 
mean of two terms representing a two-stage process, in which a voter in one track converts 
an individual of the opposite affinity class into a third-party voter, who then influences an 
individual in the first track to join the voting class of the original voter, thereby completing 
the cycle. Because this cross- affinity cycle has two stages, the appropriate measure of its 
efficiency is a geometric mean of the two individual stages. In the case that there is no cross- 
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affinity influence (a = 0), the tracks decouple completely at this stage, and (1271) reduces to 
R[ = max(r HH ,r LL ). 

The PFE is locally stable when R[ < 1, and unstable when R[ > 1. In other words, each 
third-party voter introduced into a population that includes high-affinity and low-affinity 
individuals must influence, on average, more than one person to vote for the third party 
during his/her voting lifetime, in order for the third-party voter classes to persist, with the 
average defined by R[. 

5.2 The Member-Free Equilibrium (MFE) and R' 2 

The two-track model also has a second threshold parameter R' 2 . Analogous to i?2 from the 
one-track model, we define R' 2 as the average number of third-party voters (V# and/or Vl) 
a member can convert into M if introduced into a population of them. Since R' 2 is primarily 
concerned with the transition from third-party voting to membership, we conveniently regard 
M as the only infectious class, in order to apply the next-generation operator method to 
determine this threshold. We then calculate, from ([5]) (or (|25|) ): 



since this is scalar we seek simply the positive part divided by the term that is subtracted 
(departing M): 



where v* H2 and v* L2 are the equilibrium values at the MFE. This expression is analogous to 
that for the one-track model, 




R 2 = ~ ( V H2 + v m) 1 
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Solving explicitly for the MFE of the two-track model is complicated, but we can show 
enough to suggest that, as expected, it is unique and exists only for R[ > 1. Since we have 
m* = 0, any MFE must be an equilibrium of the subsystem 



dt 
dv L 



=Ph(vh + crv L )(p - v H ) ~ 4>h(p ~ Vh + cr(l - p - v L ))v H - {n + e H )v H , (28) 
--(3l{ov h +v L )(l-p- v L ) - <f> L (l - p - v L + a(p - v H ))v L - (/jl + e L )v L , (29) 



dt 

since we can now rewrite h = p — vh and I = 1 — p — v^. The two resulting equilibrium 
conditions can be simplified to a single equation of degree 4, which admits up to 4 solutions. 
One of these is the PFE, which can be factored out to leave a cubic equation for the MFE. 
It can be shown that the constant term in this cubic equation is zero precisely when R[ = 1, 
so that the number of positive solutions changes by one when R[ crosses 1. 

In the special case that a = (no cross-affinity influence), the system (T2"B"|) - (l2"9"j) decouples, 
yielding equilibria Et(0, 0), E 2(1 {vh, 0), E 2b {0, v* L ), and E 2c (v* H , v* L ) , where 

v* H = p — ^ — €h and v* L — (1 — p) — 



Ph - <Ph " Pl~ 4>l 

are meaningful only if positive, i.e., if thh > 1 and r^i > 1, respectively. A straightforward 
calculation of the Jacobian matrix shows that, within this subsystem, 

• if R[ < 1 (thh < 1 and tll < 1), the only equilibrium is E\, which is locally asymp- 
totically stable (LAS); 

• if r HH > 1 and r LL < 1, Ex is unstable but E 2a is LAS; 

• if r HH < 1 and r LL > 1, E\ is unstable but E 2b is LAS; 

• if r HH > 1 and r LL > 1, E\, E 2a and £"26 are all unstable but E 2c is LAS. 
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Regardless of the value of a, it is straightforward to show that all solutions of (I28l) - (l29l) 
which begin within [0, p] x [0, 1 —p] remain within those bounds, by observing that dvu /dt < 
when vh = P and < Vi < 1 — p, and that dv^jdt < when vl — 1 — p and < v h < P- 
One can also exclude limit cycles from solutions of ff28l) - fr29l) under the assumptions that 
Ph > 4>h and (3l > 4>l, via the usual application of Dulac's Criterion: 

d ( 1 dv H \ 8 (j_dVL } <Q _ 



dvn \VhVl dt J 8vl \vhVl dt 

We can therefore apply the Poincare-Bendixson Theorem to conclude that the equilibria of 
this system with a = identified above as locally stable, are in fact globally stable. 

We can also differentiate the equilibrium conditions for (12g|) - (T29|) implicitly by a to see 
what happens as a increases from zero: for instance, 



dv i 



da 



= -PlvU 1 -P-v*l) + <Plv* l {p - v* H ) 
(p L -4> L )(l-p-2v* L )-(n + e L ) 



so 



dvi 



da 



-/3 L v* H {l-p) 



E 2a ,a=o (f3 L -4> L )(l-p)-(ii + e L )' 
Thus when run > 1 and > 1, the numerator is negative and the denominator is positive, 
so that Eia exits the state space as a increases from zero. A similar calculation holds for 

E2b- 

Since (!28l) - (f29j) is a subsystem of (pT|) - (!25|) (in which m{t) = 0), stability in the subsystem 
does not imply stability in the full system, but instability in the subsystem does imply 
instability in the full system. Thus the full system (T2TT) — (]23T) has at most one stable MFE 
when a = (and, by continuity using the above result of implicit differentiation, for a 
sufficiently small), that stability depending upon the additional dimension (m) not present 
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in fl28l) - (!29j) . as measured by R' 2 . While R 2 is not expressed explicitly, given the implicitness 
of v* H2 and v* L2) the observed uniqueness of the MFE allows us to draw conclusions about 
the two-track model numerically. 

5.3 Survival Equilibria 

Although the equilibrium conditions for (T2TT)- (l25l) are too complicated to solve outright 
(apart from the fact that v* H + v* L = for any survival equilibrium, from ( 1251) ). we can 
verify numerically not only their existence as expected when it^ > 1, but also the existence, 
under certain conditions, of a backward bifurcation at R 2 = 1 just as observed for the 
one-track model. 

Figure [6] shows a situation analogous to that depicted in Figure H] for the one-track 
model, in which survival equilibria may exist even below the PFE/MFE threshold (R2 = 
or R' 2 = 0); in Figure [7] the critical (minimum) value of R' 2 is between and 1. Both figures, 
however, also demonstrate the existence of multiple locally stable survival equilibria near 
R 2 = 1 (see the close-up in Figure [6]), meaning that for some parameter values there are 
two different levels at which the party may stabilize, in addition to the stable (for R' 2 < 1) 
extinction equilibria PFE/MFE. In these situations, the initial number of party members 
plays a huge role in determining whether the party surges to major growth, languishes, or 
dies out entirely. 

Although the parameter values used to create these figures are idealized for illustration 
purposes, the value p = 0.1 reflects an estimated 10% of the voting population having high 
affinity for the Green Party's agenda, and other values reflect the distinctions in peer-driven 
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Figure 6: A bifurcation diagram for the two-track 
model illustrating a backward bifurcation at R' 2 = 1 
(see closer view at right), multiple saddle-node bi- 
furcations, and survival equilibria extending below 
R[ = 1 (R' 2 = 0). Parameter values used are p — 0.1, 
Pl = 0.2, e H = 0.25, e L = 0.5, <p H = <f) L = 0.1875, 
a = 1.25, a = 0.8, 7 = 0.6, fj, = 0.05; f3 H varies. All 
rates are in years" 1 (p, a, a are dimensionless) . 



Figure 7: A bifurcation diagram for 
the two-track model illustrating a 
backward bifurcation at R f 2 = 1. 
Parameter values used are p = 0.1, 
Ph = 10, p L = 0.2, e H = 0.25, 
e L = 0.5, <p H = <Pl = 0.1875, 
a = 1.25, a = 0.8, fj, = 0.05; 7 
varies. All rates are in years" 1 (p, 
a, a are dimensionless). 
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behavior described by the model for the two affinity classes. In addition, similar curves 
can be obtained using a wide range of values for model parameters. Unfortunately, the 
complexity of the two-track model prevents an explicit calculation of a quantity analogous 
to i? 3 which measures the ability of the member class to recruit "susceptible" voters from 
both affinity classes to become third-party voters. The same interpretation applies, however: 
it is the work of party members in recruiting voters that enables a party to persist when 
third-party voters' influence is too weak. 

6 Conclusion 

The models described in this study investigate in mathematical terms the consequences of 
our assumptions about the factors driving the growth and persistence of a third political 
party arising through grassroots efforts. In particular, we assume a hierarchical structure in 
which party members (activists) play a different, more extensive role in party survival than 
party voters. We also assume that the dominant influences are primary contacts among 
third-party voters and members and the general public, manifested in our models as nonlin- 
ear terms involving the sizes of the two groups making contact. These nonlinearities govern 
the behavior of the models — that is, the fate of the party under study — through threshold 
quantities that describe the system's tipping points. Our results should be taken as impli- 
cations of the assumptions that such primary contacts, and not other factors (apart from 
affinity as defined for our two-track model), drive individuals' decisions whether to support 
third political parties at any level. 

The primary result of our analysis is that our models identify, and provide a way to 
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measure, the three factors that determine the party's survival. Each of these factors is a re- 
productive number that describes the ability of a given class in the party structure to recruit 
others. R\ and R[ give the average number of unaffiliated voters recruited per third-party 
voter, for the simplified and general model, respectively. This first threshold quantity mea- 
sures the voting class's ability to replace or sustain itself. R2 and R' 2 , meanwhile, describes 
the efficiency with which third-party members convince third-party voters to become mem- 
bers. This average number of new members recruited per existing member presupposes the 
success of the first stage of recruitment: existing third-party voters recruiting new ones; that 
is, Ri > 1, or else R 2 < 0, which is meaningless except to indicate the failure of a recruitment 
structure in which party members play no role in recruiting new voters. Finally, -R3 quantifies 
the ability of party members to recruit new voters directly from the unaffiliated (with this 
party) public. The growth of the Green Party in states like Maine and Maryland, and its 
recent decline in states like California and Oregon, illustrate the effect of these tipping-point 
thresholds (cf. Gladwell, 2000). 

The model's prediction of survival states in scenarios where the primary, hierarchical 
recruitment structure is not strong enough (R 2 < 1, and even i? x < 1) to allow a party to 
arise — that is, political trends are not favorable and the resulting influence of peer pressure 
to vote for and support a given party is weak — explains the persistence of established third 
parties during periods of adverse conditions, when political winds blow against them. In 
particular, when party activists (members) are sufficiently capable of finding and recruiting 
new voters directly (as measured by R 3 > 1), a large enough core of committed party 
members can ensure the party's survival. In cases where the political environment was 
favorable for a time (R 2 > 1), this minimum core size is typically reached quickly. This 
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reaching across traditional hierarchical structures (rather than party members interacting 
primarily with those who already vote for the party) provides a robustness to the phenomenon 
that manifests mathematically in the backward bifurcations illustrated in earlier sections. 
Backward bifurcations also underscore the importance of initial conditions (having enough 
initial party members) in enabling that robustness. This scenario is typified by the case 
study of the Green Party in Pennsylvania. 

Our general model classifies the general voting population by affinity to the ideas and 
goals of a given third party. The form of the expression for R[ illustrates how a party's 
ability to take hold in even a small subset of the population (the high-affinity track) affects 
the party's survival in the population as a whole: since R[ is greater than either of the within- 
track voting replacement numbers (analogous to Ri), a successful enough recruitment within 
the high-affinity track can maintain the party. Furthermore, the stratification into two tracks 
creates the potential, as illustrated in Figures E] and U\ for multiple stable survival states, 
which facilitate party growth since the additional states require fewer initial members than 
the original one. Our models could easily be extended to a stratification with intermediate 
levels of affinity, each of which would make its own differential (i?i-like) contribution to 
determining the overall recruitment potential of a given third party (analogous to R[), in 
addition to increasing the diversity of possible survival states. 

Even though we used available data to estimate parameters, we cannot claim to have mea- 
sured the strength of the various influences directly. Our models provide qualitative measures 
for the efficiency of parties' recruitment strategies, identifying which factors interact, and 
how. The accuracy of these measures hinges, of course, on the underlying assumptions ar- 
ticulated in Section 12.11 Since it is often difficult in practice to quantify the strength of 

43 



another's opinions and arguments in influencing one's opinion, our models can also be inter- 
preted as an illustration of how individual interactions within a population combine to exert 
a single collective influence observable only at the population level. Finally, it should be 
noted that the deterministic nature of our models ignores the stochastic aspect of individual 
interactions (such as an exceptionally charismatic individual): individual variability is a crit- 
ical factor when considering very small groups. Here we have described political behavior in 
terms of population averages, but the first stages of any movement are entirely dependent on 
the particular personalities involved. Therefore our models should be considered as picking 
up at the point where a grassroots movement has become sufficiently organized to become 
a political force. 
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A Equilibrium Analysis for the One- Track Model 



A.l Proof of Proposition 1: Existence of E% and £4 

We begin with the survival equilibrium condition f(m*) = m* 2 + BM* + C = 0, where (from 

b = ( a _ d _ a\ + a c = -- A - ^ - ^ +e " 

\ a/3 / 7 \ 7/ a/3' 7 a/3 \ 7/ a/3 

We rewrite these expressions by defining the following terms: 

(3 — <p a af3 7 

Q = — 5-, r= — — , x = — — , y = -. 
aft jj + e jj + e jj 

By assumption, (3 > </>, 7 > //, and a > 1, so that < g < 1, < r < 1, and y > 1. In these 
terms, i?i = gx, # 2 = y (l - -^), and 

















gx / 



In order to be meaningful, solutions must fall within the interval (0, 1 — v*) (since m* > 
1 — v* will make s* = 1 — v * — m* < 0). We calculate 

(1 \ r 1 — r 
1-- =- + >0 
yj x xy 

and 

f'(l - v*) = 2(1 - v*) + B = 1 - — — + - > (since 1 - q < 1< y). 

y x 

These two inequalities imply that any roots of / lie to the left of 1 — v*. Next note that 
/(0) = C = — p-(i?2 — 1)- Then when R 2 > 1, f(0) < 0, so there is exactly one solution m\ 
in (0, 1 — v*). When R 2 < 1, /(0) > 0, so the number of solutions in (0, 1 — t> *) is even. In 
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this case, both solutions m* ± are in (0, 1 — v*) precisely when B 2 — AC > (the solutions are 
real) and B < (the parabola's vertex m = —B/2 lies to the right of 0). 

The condition B < can be shown equivalent to x > r and y > y B = (1 + O)-^.- The 
condition C > (i.e., R 2 < 1) can be shown equivalent to qx < 1 or y < y c = ^frj- The 
condition B 2 — AC > becomes 



Solving the quadratic inequality in y, this means y must not be between the two positive 



^ = { K 1 " 9)0* " + 2(1 - gr)] ± 2^(1 - qr) [{I - q){x - r) + {I - qr)]] . 

Thus in order to have two survival equilibria, we must have x > r, ys < y < yc, and either 
V < y- or y > y + . 



In order to simplify these criteria, we compare the threshold values for y. We find that 
ij- < ys < y+ is equivalent (after substitution) to 



- y/(l-qr) [(1 - q)(x - r) + (1 - qr)} < qx - 1< y/(l-qr) [(1 - q) (x - r) + (1 - qr)], 

(30) 




and multiplying by x 2 y 2 we get 



(x - r) V - 2x [(1 -q)(x-r) + 2(1 - gr)] y + (1 - g) V > 0. 



roots 



or simply |gx — 1| < a/(1 — gr) [(1 — q){x — r) + (1 — gr)]. Squaring both sides and simpli- 



fying leads to the compound inequality 



r < x < x + = 



1 + g(l -r) 

g^ 
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In this region both E 3 and E 4 exist. 




Figure 8: Four threshold curves illustrating the survival equilibrium conditions 

The lower bound r corresponds to the vertical asymptote in x shared by y± and y#. Since 
q < 1, - < 4j < x + , so at the upper bound qx — 1 > 0, and it is the second inequality in 



fl30p that is violated. That is, when x = x + , ys = y+- Thus it is always true that y- < i/B, 
so the criterion y > y B allows us to discard the condition y < y^ in favor of y > y + . 

If we solve the inequality ys < yc f° r x (when qx > 1), we find again the condition 
x < x + , indicating that the graphs of ys, yc, and y + all cross at x = x + . Further similar 
computation can show that y + < yc except at x = x + , where they are tangent. We thus 
require r < x < x + and y + < y < yc (the latter inequality only for qx > 1). A graph 
illustrating all four curves is given in Figure [HI 

To put the survival equilibrium conditions back in terms of the original parameters, we 
see that r < x < x+ becomes 



H a/3 1 + 
< < - 



/i + e /i + e 



2 ' 



51 



each of the inequalities can be solved for a/3/fi to make 



ap ( (3-<p(3-<p-e 
— > max 1 , 



Note that the last expression in the inequality above is the product of Ri and another 
fraction whose value exceeds 1 precisely when Ri does. Thus if R\ < 1, the condition is 
simply 2e > 1, while if R\ > 1 the condition is ^ > Ri((3 — — e)//i. It is not so simple 
to rewrite y > y + ; however, we can return to the original conditions B < 0, B 2 — AC > 0, 
which can be rewritten more simply, either as R% a > 1, R^ > 1 as in the statement of 
Proposition 1, or as: 



1 / 0-(j>\l 1 



IV - \C > *=>,/ A: 




a (3 



l - 1£) I + ft(f/7) < 71, 



(31) 
(32) 



where 



In order to rewrite the condition 5 2 — 4C > in the form (|32|) . we substitute S and C 
from (TT51) and get: 



a(3 



1 + 



a/3 / 7 



4 



7 



/i + e (3-4> ( // 



a/3 a/3 



7 



> 0. 



Expansion, summing like terms, and completing the square yields 



ap 



af3 J 7 



-i 2 



>4ii i + i 

a/3 V 7 



Next we take the square root of both sides, 



a/3 



- 1 



(3 — </>\ // 



a/3 / 7 



+ 1 



> 24 



a/3 
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Since we require /x < 7 and (3 > 0, then 



jlA^V + ^jl + ^ + A^, . >0 . 

ap V ap I 7 ap ap 7 V 7 



so we can drop the absolute value bars. Finally, rearranging, we get 



af3 J 7 /i y a/3 /i \ 7/ a/3 

In the case that e = 0, this inequality can be further simplified by factoring the right-hand 
side as a perfect square, taking the square root of both sides, and using the fact that fl3Tl) 
implies l/a/3 < to get 



a/3 + v( 1 a/3 ) 7 " \V 




We can apply the same technique to obtain (|32|) . 



A. 2 Stability Analysis for E% 



From the end of Section 13.31 it remains to show that tr Ji(E 3 ) < when R 2 > 2. This 
condition on the trace is equivalent to 

(/?-</>) + 

m * > T = V 

+ (P - <j> + a(3 + 7) 

In terms of q, r, x and y as defined in the previous section, 

l^-'-J), ^(i±i±S-i, =-£f 1 -±_iy 

(l + g + ^-j 27 y\qxy) 
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Now the stability condition is 



m , 



-B + VB 2 - AC 



> L 



VB 2 - AC > B + 2L 



B + 2L<0, or B + 2L > and B 2 - AC > B 2 + ABL + AL 2 



B + 2L>0 —C > BL + L 
B 2L > => -fl | >-(!- — --) -L + L 2 



V \ qx yj y V qx y 

B + 2L>0 0> — % -L + L 2 

V 



Here we multiply by — (l + q + /q, expand, and simplify to get 



< 



l + o r 
+ - 

y x 



i-l- 2 -) (i + l±*r + H? 

qx yj \ x y 



which is true since 



qx y t 

This completes the verification that E$ is LAS when it exists. 



R 2 > 2 I 1 - — - - I > 0. 



B R[ for the Two- Track Model 

We calculate the reproductive number R[ of the two-track model using the next-generation 
operator method where R[ is analogous to R\ of the one-track model. The party-free equi- 
librium E 1 of (EH) ([25D is (p, (1 -p), 0,0,0). Differentiating (12g]) -(f25l) with respect to the 
"infective" variables vh, vl, and m, and substituting the PFE values yields the following 
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"mini-Jacobian" matrix: 
/ 



A 



(Ph ~ <Ph)p ~ w H p H ap ap H p 

Pl<t(1-p) (Pl - <Pl)(1 - p) - w l ap L {l-p) 



\ 











where % = /x + e# + 0i/cr(l — p) and iol = /i + + 4>l&P- We next rewrite A 
where the entries of M are nonnegative and D is a diagonal matrix: 



= M 



( 



M = 



(Ph ~ 4>h)p Ph&P oipuP 

Pl<t(\-p) (P l - <p L )(l - p) aP L (l-p) 




and 



D = 



/i + e H + 0ho"(1 - p) 

/i + e L + (f)L<Jp 



V 



o 



o 



Now R[ is the dominant (largest) eigenvalue of 



MD V = 



{I3h-4>h)p 



fj,+€H+4>H<r(l-p) n+e L +<p L ap 



OtPHp/H 



\ 



o o 

The three eigenvalues are and 



thl r LL aP H p/ft 




thh + r LL ± a/ {r HH - r LL ) 2 + Ar HL r LH 



R[ takes the positive square root in the latter expression. By inspection we can see that 
the extreme cases p = and p — 1, R[ simplifies to R\ for the one-track model. 
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Since the second term inside the radical is positive, we have that 



Thh + r LL + a/ (r HH - r LL ) 2 



2 V'hh + r LL + \r HH - r LL \] = max(r HH ,r LL ). 



Since, for positive numbers a and b, y/cT+b < y/a + Vb, we also have that 



K<- 2 



thh + r LL + a/ {r HH - r LL ) 2 + 2y/r HL r LH j = m&x(r HH , r LL ) + ^r^Jhu- 
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Tables 
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Table 1: Compartments and parameters of the two-track model 



Table of Compartments and Parameters 


H 


high affinity susceptibles (i.e., voters highly susceptible to third party ideology) 


L 


low affinity susceptibles (i.e., voters barely susceptible to third party ideology) 


V H 


third party voting individuals deriving from H 


Vl 


third party voter individuals deriving from S 


M 


third party members (i.e., party officials, donors, volunteers) 


V 


proportion of the voting population N entering H 


Ph 


peer-driven recruitment rate of H into Vh by individuals in Vh, Vl and M 


£h 


linear recruitment rate of Vh back into H via secondary contacts 
(i.e., media and campaigning from opposing parties) 


<f>H 


recruitment rate of Vh into H by direct contact with individuals 
in the opposition classes (i.e., individuals in H and L) 


Pl 


peer driven recruitment rate of L into Vl by individuals in Vl, Vh, and M 




linear recruitment rate of Vl back into L via secondary contacts 
(i.e., media and campaigning from opposing parties) 


<\>L 


recruitment rate of Vl into L by direct contact with individuals 
in the opposition classes (i.e., individuals in H and L) 


a 


factor by which the influence of party members M in recruiting voters 
in H and L into Vh and Vl exceeds that of voters in Vh and Vl 


a 


factor by which the influence of individuals upon members of a different 
affinity class is reduced in encouraging or discouraging third party voting 


7 


recruitment rate of Vh and Vl into M by individuals in M 


fi 


rate at which individuals enter or leave the voting system 
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Table 2: Parameters of the one-track model 





peer driven recruitment rate of S into V by third party voters and members 


6 


1 t;L,l Ul Llildllj id It; Ol V 7 Udtrv 1111U O Vld oCL^UllLicll j L.U111CIL, to 

ii p mpni^ pun pa mn^ i oti in p* from nnnn^ino' nprtip^i 


6 


rpprn itmpnt, rat.p of V into "hv Hirppf, rnntapf, with siispptvHIYIps 


a 


factor by which the recruitment rate of S into V by third party 
members exceeds the recruitment rate by individuals in V 


7 


recruitment rate of V into M by third party members 


n 


rate at which individuals enter or leave the voting system 
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Table 3: Equilibria of one-track model 



n ,a mil hyi 1 1 m 


J__/-A.1oLj"11L/v; V. Will t . 


Stahilitv C!nnrl 


rarty-rree rj\ — ^i,u, uj 


always exists 


ill < 1 


Member-Free E 2 = (±, 1 - ±, 0) 


i?i > i 


R 2 < 1 


Survival E 3 = (1 - ^ - m% m^) 


(i) i? 2 > 1, or 

(ii) i? 2 < 1, R 3 > 1 


always stable 
when it exists 


Survival £! 4 = (1 - £ - ml, £ ml) 


i? 2 < 1, R 3 > 1 


always unstable 



Table 4: Regions of Equilibrium Stability 





E x 


Ei 


E 3 


E 4 


I 


stable 


unstable 


does not exist 


does not exist 


II 


unstable 


stable 


does not exist 


does not exist 


III 


unstable 


unstable 


stable 


does not exist 


IV 


unstable 


stable 


stable 


unstable 


V 


stable 


does not exist 


stable 


unstable 
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Table 5: Estimated model parameters, for the seven U.S. states (and district) used in the 
case study. States are listed in decreasing order of proportional membership M to facilitate 
comparison with Figure [5j Initial conditions v(0) marked with asterisks *, as well as all 
model parameters (except //), were estimated by data fitting as described in the main text. 
All rates are given in units of 1/days; values below 0.01 are given in scientific notation. 





Initial 


Final 


Initial Conditions 












State 


time 


time 


(5(0),V(0),M(0)) 




a(3 


€ 




7 


ME 


Jun 1998 


Nov 2006 


(91.12, 8.58, 


0.31) 


5.27E-4 


0.5022 


1.00E- 


-5 


1.34E-4 


DC 


Jun 2003 


Sep 2004 


(87.54, 8.61, 


3.85) 


0.0181 


1.9911 


1.02 




6.72E-4 




Jul 2005 


Aug 2006 


(79.52, 16.43*, 


4.05) 


1.06E-3 


0.200 


0.468 




0.0115 




Jan 2007 


Nov 2008 


(89.04, 6.63*, 


4.33) 


1.66E-3 


1.15 


2.83 




1.02E-7 


CA 


Feb 1999 


Feb 2004 


(94.87, 3.09, 


2.03) 


4.94E-4 


0.7188 


1.69E- 


-5 


3.73E-4 




Oct 2004 


May 2008 


(91.87, 4.80*, 


3.32) 


0.0105 


0.2193 


0.162 




1.94E-6 


OR 


Jan 2001 


Oct 2004 


(88.45, 10.26, 


1.29) 


1.56E-4 


1.45 


2.99E- 


-5 


1.29E-4 




Nov 2004 


Mar 2009 


(96.06, 1.70*, 


2.23) 


0.0524 


0.781 


2.93 




1.00E-6 


NY 


Apr 1999 


Mar 2004 


(99.27, 0.73, 


0.01) 


0.0119 


0.352 


4.00E- 


-3 


4.33E-4 




Nov 2004 


Mar 2008 


(98.65, 0.16*, 


1.19) 


0.0567 


0.283 


3.54 




1.00E-6 


MD 


Aug 2000 


Mar 2009 


(98.17, 1.32, 


0.01) 


0.0103 


1.03 


2.40E- 


-3 


1.93E-4 


PA 


Apr 2001 


Nov 2006 


(97.33, 2.55, 


0.11) 


1.01E-5 


0.329 


1.37E- 


-6 


9.85E-4 
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Table 6: Reproductive numbers and stable equilibria for the one-track model, calculated 

using parameter estimates from Table [51 Reproductive numbers not given are negative. 

Equilibria are given as percentages of the target population. 

State Period R x R2 R 3 Stable equilibria (v*, m*) 

ME 1998-2006 9.44 2.62 43.5 E 3 = (34.15, 65.84) 

DC 2003-2004 0.18 - 72.4 E x = (0, 0), E 3 = (6.81, 89.30) 

2005-2006 0.0023 - 38.9 E x = (0, 0), E 3 = (0.39, 98.63) 

2007-2008 0.00059 - - E 1 = (0, 0) 

CA 1999-2004 7.88 7.11 81.3 E 3 = (12.28, 87.71) 

2004-2008 0.065 - - E x = (0, 0) 

OR 2001-2004 2.06 1.45 71.8 E 3 = (35.40, 64.59) 

2004-2009 0.018 - - E x = (0, 0) 

NY 1999-2004 2.93 6.23 53.4 E 3 = (10.58, 89.27) 

2004-2008 0.016 - - E x = (0, 0) 

MD 2000-2009 4.20 3.21 72.1 E 3 = (23.76, 76.16) 

PA 2001-2006 0.21 - 64.8 E x = (0, 0), E 3 = (4.84, 95.14) 
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